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Time-Domain  System  for  Identification  of  the  Natural  Resonant 
Frequencies  of  Aircraft  Relevant  to  Electromagnetic  Compatibility  Testing 

John  W.  Adams,  Arthur  R.  Ondrejka  and  Herbert  W.  Medley 

A method  of  measuring  the  natural  resonant  frequencies  of  a 
structure  is  described.  The  measurement  involves  irradiating  this 
structure,  in  this  case  a helicopter,  with  an  impulsive 
electromagnetic  (EM)  field  and  receiving  the  echo  reflected  from  the 
helicopter.  Resonances  are  identified  by  using  a mathematical 
algorithm  based  on  Prony's  method  to  operate  on  the  digitized 
reflected  signal.  The  measurement  system  consists  of  special  TEM 
horns,  pulse  generators , a time-domain  system,  and  Prony's  algorithm. 
The  frequency  range  covered  is  5 megahertz  to  250  megahertz;  this 
range  is  determined  by  antenna  and  circuit  characteristics . 

The  measurement  system  is  demonstrated,  and  measured  data  from 
several  different  helicopters  are  presented  in  different  forms. 
These  different  forms  are  needed  to  determine  which  of  the  resonant 
frequencies  are  real  and  which  are  false.  The  false  frequencies  are 
byproducts  of  Prony's  algorithm. 

Key  words:  avalanche  step  generator;  natural  resonant  frequencies; 

Prony's  method;  TEM  horn  antennas;  time-domain  measurement  system. 


1 . 0 Background 


If  a structure  is  exposed  to  electromagnetic  fields,  resonances  enhance 
the  fields  by  a factor  of  Q and  often  cause  electromagnetic  compatibility 
(EMC)  problems.  Structures  such  as  helicopters  are  built  of  highly 
conductive  materials  such  as  aluminum,  which  typically  yield  Qs  of  10  to 
200.  The  quality  factor,  Q,  is  the  ratio  of  stored  to  dissipated  energy. 
If  the  shape  of  a structure  is  complicated,  as  is  the  case  with  helicopters, 
these  resonant  frequencies  cannot  be  calculated.  The  only  way  to  obtain 
these  resonances  is  by  measurements.  Fortunately,  advances  in  antenna 
design  and  solid  state  technology  can  be  used  with  existing  mathematical 
procedures  to  develop  a measurement  system  that  can  make  the  needed 
measurements . 

A technological  advance  that  will  change  some  of  the  aircraft 
character i st i cs  is  the  increasing  use  of  composite  materials  in  place  of 
aluminum.  For  lossless  structures , size  - frequency  scaling  might  be  used  for 
the  type  of  measurements  discussed  in  this  paper,  but  for  lossy  structures, 
measurements  should  be  made  on  full-size  aircraft  rather  than  scale  models, 
since  size-frequency  scaling  does  not  work  effectively  for  lossy  objects. 

Prony's  algorithm  has  been  in  existence  for  about  200  years.  Morf 
recent  applications  to  determine  poles  from  time-domain  waveforms  have  met 
with  varying  degrees  of  success  [1-5].  As  a tool  for  target  identification, 
it  has  only  limited  success,  since  the  quality  of  results  degrades  as  the 
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si gnal -to-noise  ratio  decreases.  If  the  target  is  not  stationary  with 
respect  to  a transmitting  antenna,  only  one  time  interval  of  data  may  be 
recorded.  However,  if  the  position  of  a test  object  is  fixed  relative  to 
transmitting  and  receiving  antennas,  a large  number  of  time  intervals  may  be 
recorded  and  averaged,  which  gives  a significant  improvement  in  signal-to 
noise  ratio. 


Time-domain  technology  has  been  more  easily  handled  in  theory  than  in 
application,  but  its  great  potential  can  now  be  realized  due  to  the  recent 
improvements  in  the  electronics  and  antennas  needed  for  its  application. 
Higher  speed  devices  such  as  gallium  arsenide  transistors  [6-7]  are  one  area 
where  significant  improvements  are  resulting  in  faster  and  more  powerful 
impulse  generation  capability.  Antennas  that  radiate  impulses  without 
significant  phase  distortion  are  fairly  straightf orward  to  construct,  and 
the  needed  theory  concerning  differences  between  transmitting  and  receiving 
horn  antennas  is  now  available  [8].  These  transverse  electromagnetic  (TEM) 
horn  antennas  have  been  built  and  evaluated  previously,  but  only  for  higher 
frequency  ranges.  For  EMC  testing,  the  lower  frequency  ranges  typical  of 
resonances  of  larger  structures  had  not  been  covered. 


Electromagnetic  compatibility  (EMC)  testing  for  incident  fields  should 
be  performed  at  all  frequencies  to  which  the  device  or  system  may  be 
exposed.  The  lower  portion  of  the  spectrum  where  large  structures  may  be 
resonant  is  often  not  tested.  This  is  because  common  test  facilities  have 
severe  limitations  which  preclude  their  use  in  these  lower  frequencies. 
Anechoic  rooms  have  increasingly  high  reflections  at  lower  frequencies.  TEM 
ceils  or  parallel  plate  structures  can  be  used  at  low  frequencies,  but  have 
practical  limitations  for  size  of  test  object.  Outdoor  ranges  cannot  be 
used  for  swept  frequency  testing,  since  use  of  field  strength  levels  high 
enough  for  immunity  tests  would  not  comply  with  FCC  regulations.  Even  for 
EMC  testing  under  these  constraints,  the  testing  can  be  much  more  effective 
if  the  resonant  frequencies  of  the  structure  are  known.  Finding  these 
frequencies  is  the  purpose  of  the  work  reported  here. 


2 . 0 The  Measurement  System 

The  measurement  system  used  at  Ft.  Rucker  on  several  Army  helicopters 
is  shown  in  figure  1.  The  aircraft  under  test  (AUT),  whose  resonances  are 
being  measured,  is  located  on  a ground  screen  which  is  continuous  from 
transmitting  antenna  to  receiving  antenna.  The  transmitting  antenna 
launches  an  impulsive  field  toward  the  AUT,  and  properly  polarized 
reflections  or  reradiations  are  detected  at  the  receiving  antenna.  The  two 
antennas  are  oriented  orthogonally  to  each  other  to  reduce  the  direct 
coupling  of  energy  between  them.  The  AUT  was  in  its  flight  configuration, 
with  no  tow  bars  connected  or  doors  or  panels  open,  but  in  a power-down 
condition . 


The  actual  resonant  frequencies  are  unknown,  so  that  it  is  most 
efficient  to  illuminate  a target  with  a spectrum  that  contains  as  many 
different  frequencies  as  possible.  The  impulse  waveform  has  a spectrum 
which  contains  a broad  range  of  frequencies.  The  spectrum  of  other 
waveforms  usually  favors  a narrower  band  of  frequencies. 
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Time-domain  antenna  theory  shows  that  an  antenna's  characteristic  in 
the  transmitting  mode  is  the  time  derivative  of  the  receiving  mode 
response  [8].  To  produce  a radiated  impulse  field,  it  is  necessary  to 
drive  the  transmitting  antenna  with  a step  waveform. 

The  outgoing  impulse  is  produced  by  a step  function  followed  by  a slow 
ramp  which  returns  the  voltage  to  its  original  value.  The  waveform  is 
repetitive  at  a rate  determined  by  the  clock  generator.  The  trigger  pulse 
and  dc  power  are  fed  to  the  step  generator  through  a combiner. 

The  receiving  antenna  output  can  be  viewed  on  an  oscilloscope,  which 
also  serves  as  a waveform  digitizer.  For  better  control  of  the  scope's  sweep 
delay,  the  scope  was  synchronized  externally  to  the  clock  generator.  The 
computer  was  programmed  to  record  the  digitized  data  so  it  could  be 
transmitted  to  a mainframe  computer  for  further  analysis  at  a later  time. 

2.1  Avalanche  Step  Generator 


The  step  generators  used  on  this  project  were  of  the  avalanche  type. 
After  testing  many  transistors  for  maximum  applied  voltage  before  avalanche, 
consistency  of  avalanche  voltage,  and  switching  speed,  we  used  a 2N390M  npn 
transistor  in  a single  stage  for  the  preliminary  test.  In  the  final 
configuration,  we  increased  the  number  of  stages  from  one  to  four  to 
increase  the  output  amplitude  by  about  a factor  of  four.  This  was  necessary 
due  to  the  low  amplitude  of  the  differentiated  output  pulse.  The  2N390M 
transistor  switches  slowly  enough  in  the  avalanche  mode  that  a substantial 
amount  of  the  applied  voltage  is  lost  before  the  pulse  can  be  formed.  This 
is  due  in  part  to  the  antenna's  discharging  before  the  transistor  can  be 
fully  switched  on.  Also  some  of  the  applied  voltage  is  lost  due  a voltage- 
divider  effect  when  the  pulse  passes  impedance-matching  resistors  at  the 
input  of  the  antenna.  These  resistors  are  necessary  to  prevent  multiple 
reflections  if  part  of  the  pulse  is  reflected  from  the  mouth  of  the  antenna. 
These  effects  caused  the  actual  output  pulse  to  be  only  3 to  M volts  in 
amplitude  at  the  mouth  of  the  horn  and  1 to  2 nanoseconds  in  duration  at  the 
transmitting  antenna,  and  the  received  pulse  with  this  configuration  was 
only  30  to  ^40  millivolts  in  amplitude. 

In  subsequent  testing,  we  found  an  experimental  transistor  that  showed 
signs  of  being  much  faster  in  the  avalanche  mode  at  substantially  lower 
applied  voltages  than  others  tested.  That  this  transistor  would  avalanche 
at  a lower  voltage  is  a convenience,  but  the  faster  avalanche  time  gives  a 
much  faster  rise  time  as  a step-function  generator.  This  transistor  has  a 
part  number  HXTR-3615.  Another  transistor,  the  HXTR-510M,  which  is  even 
faster,  operates  at  MO  to  100  volts  dc  with  specified  switching  speeds  of  20 
to  80  picoseconds.  In  the  single  stage  version,  this  transistor  achieved  an 
output  pulse  of  M50  picoseconds  with  an  amplitude  of  M0  volts.  The 
switching  times  are  greatly  reduced  by  shortening  lead  lengths  and  by 
selective  screening  of  components. 

The  step  generator  circuit  shown  in  figure  2 is  similar  to  the  one  that 
drives  the  transmitting  antenna,  but  can  be  modified  to  operate  without  th-..- 
antenna  attached.  The  combined  dc  high  voltage  and  a 5-volt  positive 
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trigger  pulse  are  applied  to  the  input.  The  dc  voltage  charges  the  220 
picofarad  capacitor  to  a voltage  just  a little  less  than  that  needed  to 
avalanche  the  transistor.  The  trigger  pulse  coupled  to  the  transistor  base 
through  the  1 000 -pi cofarad  capacitor  initiates  the  avalanche  with  very 
little  time  jitter.  For  test  purposes,  the  step  output  is  taken  across  the 
50 -ohm  resistor  at  the  output.  In  actual  use,  the  220-picofarad  capacitor  is 
replaced  with  the  antenna  which  has  about  the  same  capacitance,  and  the  50- 
ohm  output  resistor  is  grounded. 

During  the  tests  performed  at  Ft.  Rucker,  we  ran  into  several  problems 
due  to  the  step  generators,  grounding,  and  the  ground  screen  itself: 

1.  Without  proper  grounding  of  the  step  generators,  the  pulse  would 
radiate  along  the  transmitting  trigger  cable  and  ground  screen. 

2.  The  ground  screen  was  made  of  chicken  wire,  wired  together  and 
staked  to  the  ground.  With  wind,  temperature  and  humidity  changes,  the 
ground  screen  would  shift  relative  to  other  physical  components,  causing 
changes  to  the  shape  of  the  radiated  pulse. 

3.  To  achieve  the  relatively  high  dc  voltages  necessary  for  operation 
of  the  step  generators  and  to  obtain  a clean  supply  voltage,  we  used  ten 
67.5  volt  dry  cells  in  series.  Even  though  the  step  generators  drew  very 
little  current,  a few  microamperes,  the  temperature  and  humidity  changes 
caused  the  available  voltage  to  drop  below  the  operating  range  of  the  step 
generators  under  some  conditions. 

To  solve  these  problems  we  have  done  several  additional  tests  and 
suggest  the  following  improvements  for  any  subsequent  tests: 

1.  Work  has  been  done  on  the  component  selection  and  mechanical 
design  of  the  step  generators.  After  testing  the  various  transistors,  we 
have  chosen  two  versions  for  a monolithic  chip.  Using  these  and  selectively 
choosing  them  for  specific  switching  and  avalanche  characteristics,  we  can 
have  them  assembled  on  a substrate  by  a hybrid  IC  manufacturer.  By  using 
hybrids  assembled  in  a machined  brass  package  with  male  and  female  SMA 
connectors,  we  can  achieve  a small,  1.27  x 1.91  x 0.6M  centimeter  (0.5  x 
0.75  x 0.25  inch),  shielded  package  that  can  be  bolted  to  the  ground  screen, 
thus  solving  our  grounding  problem.  The  step  generator  can  be  optically 
triggered,  further  reducing  cable  radiation. 

2.  In  the  immediate  vicinity  of  the  antennas  we  would  construct  i 
rigid  ground  screen  from  three  sections  of  0.91  by  1.22  meter,  0.16- 
centimeter  thick  printed  circuit  board  material  with  2 ounce  copper  per 
square  foot  (28.35  grams  per  0.0929  square  meters)  on  one  side.  We  would 
join  these  with  10.16  centimeter  wide,  adhesive  copper  foil  tape.  This 
would  give  a 0.91  by  3.66  meter  rigid  mounting  surface  for  the  step 
generators  and  antennas.  This  would  stop  the  movement  of  the  antenna 
relative  to  the  ground  screen. 

3.  Although  the  new  step  generator  design  allows  much  lower  operating 
voltages,  we  feel  that  a well  regulated  and  filtered  power  supply  would  be 
the  best  choice  to  power  the  step  generators.  We  have  located  a power 


supply  that  delivers  600  volts  dc  at  1.5  amperes  and  has  remote  current  and 
voltage  sense  features.  It  is  also  IEEE  488  bus  controllable.  This  source 
would  also  eliminate  downward  dc  voltage  drift  typical  of  batteries  being 
drained  by  use. 

All  of  the  above  equipment  works  well  with  the  half  horn  antennas  which 
are  suitable  for  radiating  pulses  with  vertically  polarized  E fields,  but 
the  full-horn  antennas  needed  for  radiating  horizontally  polarized  E fields 
suffered  signal  loss  due  to  ground  loading  when  they  were  in  close  proximity 
to  the  ground  screen.  To  reduce  this  effect,  tripod  supported  stands  to 
raise  the  antennas  3.7  to  4.6  meters  above  the  ground  are  needed.  These 
stands  would  be  adjustable  from  the  ground  in  both  the  horizontal  and 
vertical  planes  so  the  antennas  could  be  focused  on  the  test  object. 

Tests  have  been  performed  with  these  new  step  generator  designs  on  both 
the  indoor  and  outdoor  ground  screens  at  NBS  and  the  results  are  very 
promising . 

2.2  TEM  Horn  Antennas 


The  antennas  used  for  these  measurements  have  to  transmit  and  receive 
impulsive  waveforms  with  minimum  distortion.  They  were  developed  from  a 
basic  TEM  horn  design  described  by  Lawton  and  Ondrejka  [8].  In  the 
receiving  mode  they  have  a usable  bandwidth  from  5 megahertz  to  1 gigahertz 
over  which  the  response  is  flat  within  ±6  decibels. 

The  antenna  is  constructed  in  two  forms,  a half  horn  that  is  mounted 
directly  onto  a ground  plane,  and  a balanced  horn  that  can  be  used 
independently  of  a ground  plane.  The  half  horns  are  useful  when  operating 
with  vertically  polarized  fields;  the  balanced  horns  can  be  used  with 
arbitrarily  linearly  polarized  EM  fields.  Most  of  the  measurements  at  Ft. 
Rucker  were  made  using  the  half-horn  configuration. 

2.2.1  Receiving  Antennas 


Figure  3 shows  the  construction  of  the  balanced  receiving  antenna.  The 
two  triangular  conductors  form  a balanced  200-ohm  transmission  line,  a 
simple  modification  of  the  parallel  plate  transmission  line.  The  spacing 
between  the  conductors  tapers  to  a point  where  they  connect  to  the  200 -ohm 
input  of  a broadband  pulse  transformer.  This  terminates  the  antenna  in  its 
characteristic  impedance,  which  minimizes  reflections.  The  output  of  the 
transformer  is  unbalanced  50  ohms,  which  matches  the  coaxial  cable  we  used. 
The  upper  frequency  limit  is  about  1.2  gigahertz  due  to  the  transformer. 
The  lower  frequency  limit  is  about  5 megahertz  due  to  the  antenna. 

When  the  antenna  is  mounted  on  the  ground  plane,  only  half  of  this 
structure  is  needed.  The  second  half  appears  as  an  image  in  the  ground 
plane.  Therefore,  we  refer  to  this  as  the  half-horn  antenna.  For 
convenience,  the  antenna  is  constructed  to  have  a 50  ohm  impedance  so  that, 
the  antenna  can  directly  feed  a 50-ohm  measurement  system.  The  upper 
frequency  limit  of  the  antenna  is  3*0  gigahertz,  while  the  lower  frequency 
limit  is  5 megahertz. 
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Other  antennas  are  relatively  broadband  but  often  do  not  have  good 
phase  response.  With  excessive  phase  shift,  an  antenna  disperses  a pulse  by 
delaying  different  frequency  components  by  differing  amounts.  This  phase 
distortion  makes  it  difficult  to  perform  radiated  time-domain  measurements 
and  obtain  meaningful  results.  The  NBS  antennas  have  a linear  phase 
response  over  the  passband,  making  possible  the  acquisition  of  a high 
fidelity  waveforms. 

2.2.2  Transmitting  Antennas 


The  transmitting  antennas  are  the  same  as  the  receiving  antennas  except 
that  a transformer  was  not  used.  Instead,  the  antenna  elements  were 
connected  directly  to  the  step  generator  described  in  a previous  section. 
The  antenna  itself  is  the  pulse-forming  capacitor  for  the  step  generator  and 
is  charged  to  the  high  dc  voltage  needed  for  breakdown. 

Figure  shows  the  half-horn  antenna  including  a step  generator.  When 
the  step  waveform  propagates  to  the  aperture,  the  antenna  produces  two  types 
of  fields,  a step-like  function  in  the  near  field,  and  an  impulse  function 
in  the  far  field.  Only  the  impulse  radiates  at  a significant  distance  from 
the  antennas. 

Not  all  the  dc  energy  stored  in  the  antenna  is  radiated.  Some  is 
reflected  from  the  aperture  and  propagates  back  to  the  step  generator  which 
remains  in  a state  of  avalanche.  The  resistor  R1  in  figure  2,  located 
between  the  transistor  collector  and  the  antenna,  matches  the  antenna 
impedance  so  that  multiple  reflections  are  prevented.  Figure  5 shows  the 
difference  in  radiated  waveform  for  various  values  of  R1  . There  are  very 
few  reflections  in  the  waveform  from  the  matched  generator  compared  to  the 
unmatched  one,  and  most  of  these  are  external  to  the  antenna.  However,  the 
peak  amplitude  in  the  matched  condition  is  half  that  of  the  unmatched 
condition  since  the  matching  termination  for  the  source  requires  half  the 
voltage. 

Because  the  received  signal  is  always  contaminated  by  noise,  it  is 
desirable  to  increase  the  amplitude  of  the  transmitted  signal  as  much  as 
possible.  A modified  circuit,  shown  dashed  in  figure  2,  introduces  a diode 
shunting  R1  to  enhance  the  output  step  while  still  retaining  the  200 -ohm 
load  resistor  for  reflections.  The  amplitude  of  the  resulting  pulse  with 
the  diode  is  50  percent  larger  than  the  matched  pulse  without  the  diode.  No 
additional  reflections  are  produced  by  the  diode. 

2 . 3 Waveform  Measurement  System 


Signals  from  the  antennas  are  measured  directly  on  a digitizing 
oscilloscope.  The  input  plugin  unit  to  the  oscilloscope  is  capable  of 
sampling  waveforms  as  fast  as  3 nanoseconds,  with  an  equivalent  bandwidth  in 
excess  of  200  megahertz.  At  these  frequencies,  the  oscilloscope  has  to  be 
operated  in  the  repetitive  sampling  mode.  The  scope  is  capable  of  multiple 
acquisitions  and  can  calculate  an  average  waveform.  It  took  about  a minute 
to  acquire  an  average  of  1000  waveforms  the  way  the  scope  was  set  for  use  at 
Ft.  Rucker.  This  amount  of  averaging  enhanced  the  si gnal-to-noise  ratio  by 
about  a factor  of  30. 
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The  digitized  information  was  also  transferred  to  a personal  computer 
via  the  IEEE-488  bus  where  it  could  be  displayed  or  transformed  to  the 
frequency  domain.  After  this  optional  processing  of  the  waveform,  the 
digitized  data  file  was  stored  on  a floppy  disk  for  additional  processing 
using  more  powerful  algorithms  on  a mainframe  computer.  The  programs  needed 
in  the  personal  computer  included  a data  transfer  routine  to  establish 
communication  between  the  computer  and  the  oscilloscope,  a program  that 
could  calculate  the  fast  Fourier  transform  (FFT)  and  graph  both  the  time- 
domain  waveform  and  the  its  spectrum,  a file  conversion  routine  that  could 
change  the  data  file  from  that  which  was  generated  by  the  oscilloscope  to  a 
form  that  could  be  transferred  to  a mainframe  computer,  and  a program  which 
could  convert  data  from  the  mainframe  computer  into  a form  suitable  for 
graphic  displays. 

2. 4 Use  of  Prony's  Algorithm 


Prony's  algorithm  is  a mathematical  procedure  to  extract  frequencies 
and  corresponding  residues  from  a time-domain  waveform  which  contains  damped 
sinusoids,  as  illustrated  in  figure  6.  After  the  input  data  is  processed, 
the  output  data  is  given  as  frequencies  and  residues  in  the  s plane.  The 
frequencies  appear  as  conjugate  pole-pairs  in  the  second  and  third  quadrants 
of  the  complex  s plane;  only  the  second  quadrant  poles  are  shown  in 
subsequent  figures,  with  radians  converted  to  hertz. 

The  number  of  pole  pairs  is  usually  unknown.  Since  one  of  the  inputs 
needed  for  Prony's  algorithm,  in  addition  to  the  time-domain  waveform,  is 
the  number  of  poles,  we  had  to  experiment  with  known  waveforms  and  try 
requests  for  a different  numbers  of  poles.  We  found  no  satisfactory 
explanation  or  guidance  in  the  literature.  A heuristic  guideline  that  we 
came  up  with  works  consistently,  but  for  the  mathematical  purists,  nothing 
short  of  a rigorous  proof  is  convincing.  In  any  case,  additional 
theoretical  work  might  be  appropriate  to  evaluate  the  procedure  we  developed 
by  trial  and  error  on  known  waveforms.  This  procedure  is  explained  as 
follows . 

Repeated  runs  requesting  different  numbers  of  poles  from  Prony's 
algorithm  for  each  data  set  are  a key  to  this  procedure.  The  number  of 
poles  is  not  known  beforehand.  If  there  are  n true  poles  and  m poles  are 
requested  where  m < n,  Prony's  algorithm  breaks  down  completely  and  gives 
randomly  located  false  poles  that  are  not  even  close  to  the  true  poles.  If 
m poles  are  requested,  where  m > n,  then  n poles  are  correctly  located, 
while  m-n  poles  are  randomly  located.  If  p poles  are  requested  where  p > 
n,  and  p * m,  n poles  are  again  correctly  located,  and  p-n  poles  are 
randomly  located.  Repeating  this  request  about  four  times  seems  to  be 
sufficient  to  tell  which  poles  have  identical  locations  and  can  therefore  be 
identified  as  true  poles.  The  randomly  located  poles  only  overlap  by 
chance,  and  in  the  two-dimensional  s (complex  frequency)  plane,  this  seldom 
occurs,  certainly  not  for  all  four  requests. 

These  conditions  are  shown  in  four  figures  produced  from  a genera*  ' • 
curve  with  25  poles  from  7 to  246  MHz.  If  m or  p are  greater  than  50  (f  r 
25  complex  conjugate  pole  pairs),  the  25  frequencies  are  shown  correctly 
figures  7 and  8 by  the  repeated  poles  thqt  form  a tilted  Z.  In  figur- 
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the  solutions  from  request  of  55  poles  (x's)  and  request  of  50  poles  (V's) 
are  correct;  the  solutions  from  request  of  45  poles  (#’s)  are  close  in  some 
cases  but  substantially  in  error  other  cases;  the  solutions  from  request  of 
40  poles  is  grossly  in  error  everywhere.  Figure  10  shows  the  meaningless 
scatter  obtained  when  all  requests  were  for  less  than  50  poles. 

Although  plots  can  give  this  information,  it  is  much  clearer  to  watch  a 
buildup  of  the  poles  on  a CRT  display  in  order  to  separate  actual  poles  from 
artificial  poles.  For  example,  it  is  hard  to  distinguish  between  two,  three 
or  four  poles  printed  at  the  same  location.  For  this  we  used  a personal 
computer  with  a high  resolution  video  display. 

We  used  this  procedure  and  Prony's  algorithm  to  determine  roots  of 
these  digitized  time-domain  waveforms.  We  generated  three  different  formats 
for  each  set  of  output  data.  These  three  formats  are  (1)  a plot  of  the 
poles  in  the  second  quadrant  of  the  s plane,  (2)  a bar  plot  showing  the 
magnitude  of  the  complex  residues  of  the  poles  and  (3)  a bar  plot  of  the  Qs 
of  the  poles. 

The  pole  plot  in  the  s plane  shows  the  approximate  natural  resonant 
frequency  of  a pole  (proportional  to  the  jw  value)  and  the  damping  factor  of 
each  pole  (related  to  the  o value).  Poles  near  the  frequency  axis  indicate 
slowly  damped  sinusoids  and  will  have  high  Qs . Poles  further  to  the  left 
indicate  more  rapidly  damped  sinusoids  and  will  be  of  little  concern  since 
their  Q values  are  low. 

The  bar  plot  of  the  magnitude  of  the  complex  residues  shows  their 
amplitudes.  The  amplitude  is  related  to  the  power  of  the  frequency  source 
for  the  particular  pole  and  is  therefore  related  not  only  to  the  gain  of 
both  antennas  of  our  measurement  system  but  also  to  the  transmitting  and 
receiving  characteristics  of  the  AUT.  The  residues  may  also  be  used  several 
ways  to  assist  in  the  separation  of  true  poles  from  artificial  poles.  The 
magnitude  of  the  complex  residue  is  one  dimensional,  so  random  repetitions 
have  a better  chance  of  occurring  than  in  the  two-dimensional  s-plane  plots 
above.  The  residues  repeat  for  the  four  runs  for  the  lower  frequency  poles 
only  if  the  algorithm  has  a true  pole.  The  residues  of  the  artificial  poles 
will  not  repeat  unless  the  magnitude  is  extremely  low.  A very  low  residue 
is  typical  of  poles  extracted  from  system  noise.  These  show  a very  high  Q 
and  may  be  fairly  tightly  clustered.  The  poles  with  low-residues  appear  to 
be  repetitive  because  visual  displays  of  very  small  numbers  seem  about 
equal  . The  high  residues  typical  of  poles  in  the  lower  frequency  range  are 
not  repeated  for  all  four  runs  and  are  therefore  relatively  easy  to  use  in 
identifying  true  or  false  poles. 

The  Q of  each  pole  is  given  by  the  ratio  of  imaginary  and  real  parts  of 
the  complex  pole.  For  EMC  applications  this  is  a very  important  parameter 
since  it  gives  an  estimate  of  how  much  unperturbed  fields  may  be  enhanced  by 
the  resonance  of  the  structure. 

3.0  Measurement  Procedure  Used  at  Ft.  Rucker 

There  are  several  steps  to  obtain  the  final  answers.  First,  the  time- 
domain  system  and  transmitting  and  receiving  antennas  must  be  set  up  in  a 
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way  that  minimizes  cross  coupling  and  reflected  signals  from  unwanted 
objects,  and  that  maximizes  desired  signals  from  the  AUT.  Second,  the  time- 
domain  parameters  such  as  the  sampling  rate  must  be  selected  appropriately 
for  the  preceding  placement.  Third,  reference  signals  must  be  recorded  with 
no  test  object  present  so  that  this  signal  can  be  subtracted  from  signals 
with  test  object  present  in  order  to  remove  the  effects  of  reflections  and 
some  time-domain  system  characteristics . Fourth,  signals  must  be  obtained 
for  several  orientations  of  the  AUT  since  coupling  varies  depending  on  the 
antenna  pattern  of  the  test  object  and  the  two  measurement  system  antennas. 
Fifth,  the  filtered  and  digitized  data  must  be  recorded  on  disks.  Sixth, 
these  digitized  data  must  be  processed  using  Prony's  method  and  the  output 
data  must  be  displayed  in  several  different  forms  in  order  to  obtain  the 
correct  interpretation. 

3 . 1  System  Setup 


The  transmitting  and  receiving  antennas  were  set  up  on  suitable  ground 
screens  with  a 90°  angle  between  them,  with  both  aimed  at  a center  point  as 
shown  in  figure  1.  Typical  distances  are  shown  but  some  adjustments  were 
necessary  depending  on  the  size  of  the  AUT.  It  was  necessary  to  use  a dummy 
target  during  setup  to  obtain  large  enough  received  signal  to  make  the 
initial  setup  of  the  measurement  system. 

3 . 2  Selecting  Appropriate  Parameters  for  Time-Domain  System 


The  sampling  rate,  the  time  between  samples  and  the  number  of  samples 
to  be  used  for  averaging  were  determined  based  on  size  of  test  object, 
separation  distances,  received  signal  strength,  and  desired  frequency  range 
to  be  covered.  The  most  commonly  used  values  at  Ft,  Rucker  were  the 
following.  The  clock  frequency  was  10  kilohertz,  the  samel ing  density  was 
set  for  512  points  per  waveform,  and  the  scope  sweep  speed  to 
100  nanoseconds  per  division.  These  settings  gave  a 1.953  nanosecond  time 
interval  between  samples. 

The  clock  frequency  was  chosen  to  be  high  enough  to  allow  acquisition 
of  a data  set  in  less  than  one  minute,  but  low  enough  that  the  average 
current  does  not  degrade  the  avalanche  transistor  due  to  overheating.  The 
scope  sweep  speed  was  chosen  to  capture  all  of  the  reflections  from  the 
aircraft.  The  512  points  per  waveform  gives  very  good  time  resolution  of 
the  reflected  waveform. 

3.3  Taking  Reference  Data 


A reference  run  with  no  test  object  present  was  made  at  the  beginning 
and  end  of  each  data  recording  session.  In  addition  to  the  measurement 
system  characteristics  that  need  to  be  removed  later  by  subtracting  this 
reference  data,  there  were  many  other  external  signals  that  were  removed 
both  by  this  subtraction  and  by  averaging.  The  strongest  signals  came  from 
the  air-to- ground  and  ground-to-air  radio  communications.  These  would 
literally  cover  the  screen  of  the  oscilloscope  at  times.  There  were  some  F" 
stations  and  other  signals  of  unknown  origin. 
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3. 4  Rotating  the  Test  Object 


The  AUT  was  rotated  through  145°  increments  for  a total  of  eight 
orientations.  The  natural  resonant  frequencies  of  the  AUT  do  not  change 
with  this  rotation,  but  the  coupling  of  energy  in  and  out  of  the  AUT  would 
vary  depending  on  the  patterns  of  the  antenna  formed  by  the  AUT  and  by  our 
transmitting  and  our  receiving  antennas.  If  any  one  orientation  couples  a 
very  small  amount  of  energy  because  of  a null  in  an  antenna  pattern, 
important  information  might  have  been  lost.  The  optimum  number  of 
orientations  needed  are  not  known  due  to  the  structural  complexity  of  the 
AUT.  The  eight  rotations  we  used  seemed  to  be  adequate. 


For  each  rotation  of  the  AUT,  requests  for  four  different  numbers  of 
poles  (99,  98,  97  and  96)  were  made  of  Prony’s  algorithm.  Data  from  a 
single  rotation  is  shown  in  the  illustrations. 

3.5  Filtering  Analog  Data  Before  Digitizing 

Analog  data  must  be  filtered  before  it  is  digitized  to  prevent  aliasing 
[1].  Using  a low-pass  filter  with  sharp  roll-off  starting  at  least  an 
octave  below  Nyquist's  folding  frequency  will  prevent  signals  at  higher 
frequencies  from  being  imaged  downward  about  this  folding  frequency  and 
contaminating  the  lower-frequency  data.  Once  the  data  are  digitized,  there 
is  no  way  to  remove  these  imaged  frequencies  or  to  tell  which  are  imaged  and 
which  are  true  frequencies. 

We  did  not  pay  sufficient  attention  to  this  problem.  One  of  the  plug- 
in units  for  the  digitizing  oscilloscope  with  built-in  filtering  was  damaged 
in  shipment;  a back-up,  lower  frequency  plug-in  unit  only  had  minimal 
filtering,  and  we  did  not  have  the  proper  filters  along  to  use  as  backup. 
Therefore,  although  the  data  illustrate  the  ability  of  the  measurement 
system  to  function  as  intended,  the  data  may  possibly  be  contaminated  by 
aliased  frequency  poles.  This  aliasing  would  not  cause  frequency-shift 
errors  in  the  poles  we  did  detect,  nor  would  it  cause  true  poles  to  be 
missed.  Any  critical  use  of  these  data  should  be  verified  by  repeating  the 
measurements  and  including  sufficient  filtering. 

3. 6 Processing  and  Interpreting  the  Data 

The  digitized  data  were  then  stored  on  disk  for  later  processing  on  a 
mainframe  computer.  Some  reformatting  of  the  data  was  necessary  before  they 
were  sent  to  a mainframe  computer.  We  used  a personal  computer  for  this 
task.  We  also  used  this  computer  and  a modem  to  transfer  the  formatted  data 
to  a mainframe  computer  and  also  to  retrieve  the  processed  results.  Again, 
additional  formatting  was  necessary  on  the  personal  computer  in  order  to 
display  the  results  in  the  three  different  formats  discussed  in  section  2.H. 
We  paid  attention  to  the  adequacy  of  the  measured  signals,  the  mechanics  of 
the  processes,  and  the  adjustments  that  could  be  made  during  measurements 
subsequently.  We  overlooked  the  possible  contamination  of  data  due  to 
aliasing. 

The  key  information  to  be  obtained  from  the  processed  data  are  the 
resonant  frequencies  and  the  Q of  each  resonance.  Both  are  given 


graphically  in  the  plot  in  the  s plane.  A number  of  false  frequencies  are 
also  displayed,  but  the  printed  plots  and  CRT  displays  show  clearly  the 
correct  ones,  indicated  by  an  exact  repeat  of  the  pole  for  the  four 
different  runs  of  Prony's  algorithm.  The  exact  numbers  are  then  obtained 
from  tabulated  data.  Bar  plot  plots  of  the  magnitude  of  the  residues  offer 
additional  help  in  identifying  false  high  frequency  poles,  which  are 
indicated  either  by  a very  small  residue  or  for  larger  residues  by  lack  of 
repetition  for  the  four  runs.  The  values  of  Q are  given  in  another  bar 
plot.  The  Q plots  can  be  used  by  EMC  engineers  to  determine  an  approximate 
enhancement  of  electromagnetic  fields  by  the  resonances. 

^ . 0 Sample  Measured  Data 

Three  sets  of  sample  data  are  given.  The  first  is  a simulation  based 
on  a generated  waveform  using  known  frequencies,  damping  factors  and 
amplitudes.  These  results  are  helpful  in  learning  how  to  interpret  the 
plots.  The  second  set  of  data  is  from  the  Apache  helicopter.  The  third  set 
of  data  is  from  the  Black  Hawk  helicopter.  There  may  be  contamination  in 
both  of  the  real  sets  of  data;  the  relatively  few  clear  resonances 
identified  should  be  verified  by  additional,  independent  tests. 

Data  from  the  measurements  on  the  other  aircraft  are  summarized  in 
tables  and  figures. 

4. 1 Simulated  Data 

Figure  11  shows  a plot  of  a computer-generated,  time-domain  waveform 
which  contains  35  poles  with  frequencies  of  7,  17,  ...97,  101,  111,  ...  191, 
214,  23*1,  ...  494  megahertz.  This  set  gives  an  example  of  contamination 
caused  by  Nyquist  folding  for  frequencies  above  250  megahertz.  All  the 
poles  have  complex  residues  with  an  amplitudes  of  0.08.  Although  the  number 
of  poles  is  known  in  this  case,  we  follow  the  procedure  from  section  2.4  and 
request  a larger  number  of  poles  than  35.  Four  computer  runs  were  made 
requesting  99,  98,  97  and  96  poles.  The  repetition  of  the  location  of  the 
true  poles  is  shown  quite  clearly,  while  the  randomly  generated  poles  do  not 
repeat  and  hence  do  not  overlay.  The  digitized  waveform  has  256  points  with 
a sample  spacing  of  2 nanoseconds.  The  Nyquist  folding  frequency  for  this 
spacing  is  250  megahertz.  All  the  poles  above  250  megahertz  fold  down  and 
contaminate  the  data  below  250  megahertz.  The  pole  at  294  megahertz  is  44 
megahertz  above  250  megahertz;  it  folds  down  to  206  megahertz.  The  pole  at 
494  megahertz  is  244  megahertz  above  250  megahertz;  it  folds  down  to  6 
megahertz.  Figure  12  is  an  expanded  view  of  figure  11.  Figure  11  goes  to 
250  megahertz  while  figure  12  only  goes  to  60  megahertz.  The  residues  of 
the  folded-down  poles  are  not  changed  by  the  folding;  this  is  shown  in 
figure  13.  A bar  plot  of  Qs  is  shown  in  figure  14.  The  computer-calculated 
values  of  frequency  are  usually  correct  to  about  four  significant  figures. 
The  computer-calculated  values  of  the  residues  repeat  to  about  the  th.rc 
significant  figure  for  true  poles,  but  vary  much  mor^e  for  the  created  poles. 
Table  1 is  a computer  printout  of  the  numbers  obtained  from  these  runs. 
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Figure  15  shows  the  s-plane  plot  of  poles  to  256  megahertz,  while 
figure  16  is  an  expanded  plot  to  80  megahertz.  The  folding  frequency  for 
the  sampling  interval  of  1.953  nanoseconds  is  256  megahertz.  Resonances 
appear  at  approximately  19,  38,  59  and  76  megahertz.  The  residue  plot  in 
figure  17  also  confirms  these  as  real  poles.  Figure  18  shows  Q values  about 
60,  200,  60  and  70,  respectively.  The  printout  of  numbers  in  Table  2 shows 
the  exact  frequencies  to  be  18.6,  38.3,  59.  1 and  77.^  megahertz.  Whether 
all  these  poles  are  real  or  are  aliased  is  not  known  since  the  raw  data  were 
not  properly  filtered  before  it  was  digitized,  but  if  there  is  any  doubt 
whether  the  Apache  helicopters  might  be  expected  to  operate  in  locations 
where  external  sources  generate  these  frequencies,  additional  investigation 
is  prudent. 

M . 3 Black  Hawk  Data 

Figure  19  shows  the  s-plane  plot  of  poles  for  the  Black  Hawk  helicopter 
for  frequencies  to  256  megahertz.  Figure  20  is  the  expanded  plot  for  poles 
to  60  megahertz,  and  shows  possible  poles  at  approximately  11,  21,  28,  39 
and  55  megahertz.  These  poles  are  not  quite  as  closely  clustered  as  those 
of  the  Apache,  but  they  do  show  up  consistently  in  data  taken  with  the 
helicopter  at  other  azimuths.  Figure  21  of  the  complex  residues  casts  doubt 
as  to  whether  these  poles  are  real  in  that  the  amplitudes  do  not  repeat  as 
closely  as  those  of  the  Apache.  Figure  22  shows  the  Qs  to  be  10,  25,  23,  50 
and  50  respectively.  The  exact  frequencies  from  table  3 are  10.9,  21.3, 
28.3,  38.2  and  53.6  megahertz.  These  data  show  some  scatter  and  hence  may 
not  be  real.  The  21.3  and  28.3  megahertz  frequencies  should  be  examined  by 
additional  measurements,  since  both  are  in  frequency  bands  authorized  by  the 
FCC  for  amateur  radio  service,  and  in  addition,  the  21.3  megahertz  frequency 
is  very  close  to  the  21. M megahertz  intermediate  frequency  which  may  be  used 
in  various  avionic  systems. 

4 . M Data  From  Other  Helicopters 


The  Chinook  we  used  for  measurements  had  one  rotor  blade  removed  for 
service,  so  the  data  were  not  taken  on  a complete  aircraft. 

One  figure  of  the  expanded  pole  plot  is  given  for  each  aircraft  in 
figures  23  through  27.  Tables  M through  8 give  data  of  the  first  50  poles, 
residues,  and  Qs  for  each  of  the  aircraft  following  each  figure.  The  same 
caveat  applies  to  all  of  these  data  that  applied  to  the  earlier  data  - some 
valid  information  is  present,  but  before  it  is  used  in  any  critical 
application,  the  measurements  should  be  repeated  with  better  analog 
filtering  to  eliminate  any  possible  contamination  due  to  aliasing. 

The  processed  data  from  the  other  helicopters  and  one  fixed-wing 
aircraft  are  summarized  in  table  9 to  identify  the  resonant  frequencies. 
The  data  for  the  OH58  were  taken  with  a larger  time  interval  between 
samples,  giving  a lower  folding  frequency  (128  megahertz),  and  may  be  more 
contaminated  by  aliasing  than  the  data  for  the  other  aircraft. 
Specifically,  resonances  are  indicated  at  3 and  5 megahertz  in  figure  27, 
lower  than  is  possible  based  on  the  maximum  dimension  of  this  aircraft. 
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5.0  Summary  and  Recommendations 


A new  measurement  system  was  used  to  determine  the  natural  resonant 
frequencies  of  helicopters  by  radiating  the  aircraft  with  short  pulses, 
receiving  the  reflected  waveforms,  and  by  using  the  mathematical  transform 
known  as  Prony's  algorithm  to  extract  the  natural  resonant  frequencies  of 
these  helicopters.  Many  improvements  have  already  been  made  in  all  parts  of 
the  system;  more  are  anticipated.  One  procedural  error  was  made  in  our 
measurements.  The  analog  data  were  not  adequately  filtered  before 
digitization.  This  casts  suspicion  on  the  data  taken,  but  the  measurement 
procedure  and  system  have  been  successfully  demonstrated  at  frequencies  and 
for  applications  not  possible  previously. 

Three  additional  recommendations  are  made  for  any  subsequent  work. 
Adequate  filtering  should  be  applied.  Full-horn  antennas  should  be  used  to 
measure  horizontal  E field  polarization,  since  the  major  dimension  of  a 
helicopter  is  in  the  horizontal  plane.  The  effect  of  the  helicopter  image 
below  the  ground  screen  should  be  evaluated  using  scale  models.  This  should 
be  done  for  both  antenna  polarizations. 
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FIGURE  1 TIME  DOMRIN  MERSUREMENT  SYSTEM 
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Figure  3. 


Balanced  T£M  horn  antenna  drawings. 
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Figure  Half-horn  antenna. 
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Figure  5. 


Signal  plots  shoving  the  difference  between  matched  and 
unmatched  antenna  response. 
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Figure  6.  Sample  time-domain  waveform. 
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Figure  7.  Complex  3-plane  plot  of  second  quadrant  of  generated  data 

with  25  actual  polea  and  request  for  99,  98,  97  and  96  poles. 
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Figure  8.  Complex  3-plane  plot  of  second  quadrant  of  generated  data 

with  25  actual  pole3  and  request  for  99,  9*1,  09  and  8M  poles. 
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Figure  9.  Complex  3-plane  plot  of  3econd  quadrant  of  generated  data 

with  25  actual  poles  and  request  for  55,  50,  fl5  and  JI0  poles. 
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Figure  10.  Complex  s-plane  plot  of  second  quadrant  of  generated  data 
with  25  actual  pole3  and  request  for  >15,  ^ 3 and  >12  poles. 
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Complex  3-plane  plot  of  3econd  quadrant  of  generated  data 
po  1 **3  . 
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Figure  12.  Expanded  3-plane  plot  of  second  quadrant  of  generated  data. 
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Figure  13.  Magnitude  of  complex  residue  of  poles  for  generated  data. 
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Figure  14.  Q values  of  generated  data. 


Table  1.  Data  from  main-frame  computer  of  simulated  signal.  Table  gives 
processed  data  of  real  and  imaginary  parts  of  complex  frequency, 
real  and  imaginary  parts  of  complex  residue,  Q,  magnitude  of 
residue,  sigma  and  frequency. 
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Figure  15.  Complex  3-plane  plot  of  necond  quadrant  of  Apache  data. 
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Expanded  3-plane  plot  of  second  quadrant  of  Apache  data. 


Magnitude  of  Complex  Residues 
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Figure  17.  Magnitude  of  complex  residue  of  poles  for  Apache  data.  Units 

of  horizontal  axis  are  frequency  in  megahertz. 
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Table  2. 


Data  from  main-frame  computer  from  Apache  signal.  Table  Rives 
DrocpsspH  riat-o  , ^ . & -L c gives 


real  and  imaginary  parts  of 
residue,  sigma  and  frequency. 
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Complex  3-plane  plot  of  second  quadrant  of  Black  Hawk  data. 
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Figure  20.  Expanded  s-plane  plot  of  second  quadrant  of  Black  Hawk  data. 
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Figure  cl.  Magnitude  of  complex  residue  of  poles  for  Black  Hawk  data. 
Units  of  horizontal  axis  are  frequency  in  megahertz. 
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Figure  22.  Q values  of  Black  Hawk  data.  Units  of  horizontal  axis  are 
frequency  in  megahertz. 


Table  3-  Data  from  main-frame  computer  from  Black  Hawk  signal.  Table  gives 
processed  data  of  real  and  imaginary  parts  of  complex  frequency, 
real  and  imaginary  parts  of  complex  residue,  Q,  magnitude  of 
residue,  sigma  and  frequency. 


Complex 

Residue 

Q 

Hag 

Sigma 

Freq 

Frequency 

Res 

(Np/ys)  [MHz) 

n 6 cl 

Imag 

Real 

I mag 

-2376*00 

3 • 2662e6 

3.963e-l 

-.30907 

l 

leO 

-2.* 

3.3 

-20*8900 

3.*379e6 

2.302e-l 

-.26632 

2 

7. 0e-1 

-2.0 

3-* 

- l 8“5B00 

3-5203e6 

1 . 67le-l 

.2*12* 

2 

5. 9e-1 

-1.9 

3-5 

-188  9'i  00 

3. 5592e 6 

1 . 61 *e-l 

.25613 

2 

5 . 7e- 1 

-1.9 

3.6 

-2036500 

7 . 7 0 9 ‘ e 6 

-3-35e-l 

-.71232 

* 

leO 

-2.0 

7.7 

-1902000 

7 . 7 * 6 6e  6 

-2. 5*6-1 

- . 66255 

* 

leO 

-1.9 

7.7 

-1895000 

7 . 656Se6 

- 1 . 66e- 1 

-.701 9 

* 

1 eO 

-1 . 9 

7.9 

-1905*00 

8 . 0593- 6 

1 .8016-2 

-.7  651  1 

* 

leO 

-1.9 

8.1 

-1056*00 

1 . 067*e7 

-* . 856-2 

-.56773 

10 

leO 

-1 .0 

10.9 

-982720 

1 . 0 6 E 1 e 7 

-* . 1 6e-2 

-.5202 

1 1 

leO 

-l  .0 

10.9 

-95091 0 

1 .0692e7 

- * . 7 6e  - 2 

-.*6935 

1 1 

9 . 8e- 1 

-1 .0 

10.9 

-893*10 

1 .0909e7 

-* . 97e-2 

. *3*75 

12 

8. 6e-l 

- • 9 

10.9 

-2627900 

1 . 6602e7 

* . 660e- 1 

.12658 

6 

9. 7e-1 

-2.8 

16.6 

-25*3300 

1 . 6S9e7 

3.036e-l  - 

. 050358 

7 

6. 2e- 1 

-2.3 

1 6.9 

-2113500 

1 .7009e7 

2.5056-1  . 

0C*6i  79 

6 

5. 0e-1 

-2.1 

17.0 

-1997000 

1 . 721  7e7 

2.2026-1 

.0360*9 

Q 

✓ 

* . 5e- 1 

-2.0 

1 7 . 2_ 

-551 170 

2 . 1 26Se7 

-2 . 22e -2 

.12637 

25 

2 . 6e  - 1 

- . 9 

21.3 

-9107  60 

2. 1265e7 

-1 .966-2 

. 1*1  67 

23 

2. 9e-l 

- 0 

• «✓ 

21.3 

-51 9090 

2 . 1 295e7 

-1  .76e-2 

.11982 

26 

2 . *e- 1 

-.8 

21  .3 

-1000200 

2.1335e7 

-.3 . 1 Se~3 

.165*3 

21 

3. 3e-l 

-1 .0 

21-3 

-6000000 

2. 5756e7 

-5-216-2 

•33773 

* 

6. 8e- 1 

-6.0 

25.6 

-561  3200 

2 . 6l  5Se7 

-8. 2*e-2 

-.112 

5 

2. Se-l 

-5.6 

26.2 

-5577500 

2 . 6*8-67 

- 6 . 1 6e  - 2 

.021  * 

5 

l . 3e-i 

-5 . * 

2 6.5 

-6000000 

2.7*35s7 

-9.63e-2  - 

.005502 

5 

l . 9e-l 

-6.0 

27.* 

-12271 00 

2. 607 5e 7 

-5.87e-2  . 

0061  802 

23 

l . le-l 

-1.2 

28.1 

-121 9500 

2.8285e7 

-*.936-2  - 

. 01 6877 

23 

1 . 0e-1 

-1  .2 

26.3 

-1236200 

2 . 837367 

-*.5~e-2 

.026723 

23 

1 . 1 6 - 1 

-1.2 

28.* 

-125—00 

2 . 5* 3s 7 

-*.27e-2 

.0319 

23 

1 . le-l 

-1.3 

20.* 

-7*7200 

3 . 6002e7 

-1 .126-2  - 

.00*336 

51 

2. *e-2 

-.7 

*8.0 

-752920 

3.81167 

-1 .2*6-2  - 

.001 672 

51 

2.5e-2 

-.8 

36.1 

-761530 

3-  El  73e7 

-1 . 32e-2  - 

.001330 

*9 

2. 6e-2 

-.8 

58.2 

-8**630 

3 • 8255^7 

-1 .506-2  . 

0005878 

*5 

3.0e-2 

-.8 

36.3 

-3067*00 

* . 32“e7 

-* . 90e-2 

-.15099 

l* 

3 • 2e- 1 

-3."1 

*3.3 

-2922600 

*.565e7 

-9.356-2 

-.11273 

15 

2 . 9e - 1 

-2.9 

*3.7 

-2928700 

* .588267 

- 1 . 2 2 e - 1 

- .10*91 

15 

3 . 2e- 1 

-2.9 

*3-9 

-2965600 

* . * 1 7 1 e 7 

- 1 . 60e- 1 - 

.076927 

15 

3. 6e-l 

-3.0 

**.2 

-31*3300 

*.737Se7 

9.5576-2 

.0*80*1 

1 5 

2.le-1 

-3. 1 

*7.* 

-2993300 

* . 77  65^7 

7.3836-2  - 

.010167 

1 6 

1 . 5e-1 

-3.0 

*7.6 

-3030700 

* . 509*e7 

7. **56-2  - 

.020773 

16 

1 . 5e - 1 

-3-0 

*8.1 

-3505000 

* .66**e7 

6.3*66-2  - 

.0718*6 

l * 

1 . 9e- 1 

-3.5 

*8.6 

-1125500 

5. S5- Be  7 

1 . 699e-2  - 

.007965 

*8 

3 . 8e-2 

-1.1 

53  • 5 

1 

_* 

0 

O' 
— J 

0 

0 

5-35  6 9e  7 

1 .*736-2  . 

0072963 

50 

3-3e-2 

-1.1 

53.6 

-10*7600 

5.  36l  6 7 

1 .5*2e-2  . 

0050369 

51 

3. 2e-2 

-1.0 

5 ^ • w 

-1152500 

5. 36£9e7 

1 .7556-2  . 

0056226 

*7 

3. 6e-2 

-1.2 

53.7 

-1*1 “900 

5.978667 

-1 . 26e-2 

.01 1*99 

*2 

3. *e-2 

-1  . * 

59.8 

39 


d 66 


CL 

CL 

CL 

CO 

r- 

CO 

<T> 

CD 

cn 

0 

# 

O 

1 

s 

. ft  8 

* 

ft 

ft 

* 

o 

ft 

o 

£ 

o 

& 

w ■ - 

o o O • o o o o 

CD  CD  ■’3*  O CN  *— 

(zhiaj)  bajj 


n 

\ 

Cl 

2 


<0 

e 

tJO 

•«— t 

CO 


40 


Figure  23.  Expanded  3-plane  plot  of  second  quadrant  of  Chinook  data. 


Table  *4.  Data  from  main-frame  computer  from  Chinook  signal.  Table  gives 
processed  data  of  real  and  imaginary  parts  of  complex  frequency, 
real  and  imaginary  parts  of  complex  residue,  Q,  magnitude  of 
residue,  sigma  and  frequency. 
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Figure  2fl.  Expanded  s-plane  plot  of  second  quadrant  of  Huey  data. 


Table  5.  Data  from  main  frame  computer  from  Huey  signal.  Table  gives 
processed  data  of  real  and  imaginary  parts  of  complex  frequency, 
real  and  imaginary  parts  of  complex  residue,  Q,  magnitude  of 
residue,  sigma  and  frequency. 
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Figure  25.  Expanded  3-plane  plot  of  second  quadrant  of  Cobra  data 


Table  6.  Data  from  main-frame  computer  of  Cobra  signal.  Table  gives 
processed  data  of  real  and  imaginary  parts  of  complex  frequency, 
real  and  imaginary  parts  of  complex  residue,  Q,  magnitude  or 
residue,  sigma  and  frequency. 
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Figure  26.  Expanded  3 -plane  plot  of  second  quadrant  of  II  50  data. 


Table  7.  Data  from  main-frame  computer  of  H~58  signal.  Table  gives 
processed  data  of  real  and  imaginary  parts  of  complex  frequency, 
real  and  imaginary  parts  of  complex  residue,  Q,  magnitude  of 
residue,  sigma  and  frequency. 
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Sigma  (Np/ps) 

Figure  27.  Expanded  s-plane  plot  of  second  quadrant  of  OV-1  (fixed  wing) 


Table  8. 


Data  from  main-frame  computer  of  OV-1  (fixed  wing)  signal.  Table 
gives  processed  data  of  real  and  imaginary  parts  of  complex 
frequency,  real  and  imaginary  parts  of  complex  residue,  Q, 
magnitude  of  residue,  sigma  and  frequency. 


Complex 
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Table  9.  List  of  Aircraft,  natural  frequencies,  and  approximate  Q. 
* Aliased  frequencies  (false). 


Aircraft 

Natural 

Frequencies  (MHz) 

Q 

Chinook  CH-47 

15.1 

12 

24.2 

1 4 

28.3 

18 

32.8 

21 

36.9 

21 

ill  .2 

23 

46. 1 

25 

50. 4 

26 

Huey  UH-1 

23.5 

17 

27.4 

20 

45.6 

30 

51 .3 

38 

Cobra  AH-1 

10.3 

10 

35.2 

17 

48.2 

44 

OH-58 

3.4* 

8 

5.4* 

1 1 

17.8 

31 

22.8 

58 

26.  4 

44 

28.9 

69 

0V-1  (Fixed  wing) 

10.5 

9 

16.8 

9 

22.8 

18 

33.2 

27 

39.9 

1 7 

50.6 

17 
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